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Abstract. A fairly complete list of Toda-like integrable lattice systems, both in the contin- 
uous and discrete time, is given. For each system the Newtonian, Lagrangian and Hamiltonian 
formulations are presented, as well as the 2x2 Lax representation and r-matrix structure. The 
material is given in the "no comment" style, in particular, all proofs are omitted. 



Introduction 

This paper contains almost exclusively a large number of formulas presented in the "no 
comment "-style. I hope that they will speak for themselves and tell once more the story 
on how much contain two simple systems: the Toda lattice and the relativistic Toda 
lattice. I have only tried to collect many closely related results in the form that makes 
the interrelations quite transparent. All proofs are omitted, because they consist of a 
direct verification. 

I present a large number of integrable Newtonian systems, both in the continuous and 
discrete time, and give for all of them the equivalent Lagrangian and, for the continuous 
time systems, also Hamiltonian formulations (generalities are given in the Sect. 1), and 
also present 2x2 Lax representations. All these systems are reparametrizations of the 
Toda and the relativistic Toda lattices or their discrete time analogs (reminded in the 
Sect. 2, 3). For these systems the "big" Lax representations are well known (they include 
N x N matrices for the A-particle lattices). So, the corresponding reparametrizations 
immediately imply the "big" Lax representations for all our systems. In principle, the 
transition from the "big" Lax representation to the 2x2 one is well understood, at least 
at the formal level: they correspond to two different ways to represent a spectral problem 
connected with the second order linear difference operator. However, such transition is 
not uniquely defined, it may be performed in infinitely many different ways leading to 
different but gauge equivalent 2x2 Lax matrices. The problem is to find the gauge 
leading to the most nice Lax matrices. In particular, it is usualy required that the matrix 
L k depends only on x k , p k , the coordinates and momenta of the fcth particle. To find 
such gauge is sometimes a nontrivial task. I give here a collection of Lax representations 
having this property. 

The last Section contains several historical remarks. 

1 Newtonian, Lagrangian, and Hamiltonian 
systems in continuous and discrete time 

In what follows we consider continuous time lattice systems in the Newtonian form: 



x k = F k (x,x) 



(1.1) 



They all may be put into the Euler-Lagrange form 



d dC(x, x) dC(x, x) 
dt dx k dx k 



= 



(1.2) 



1 



We will call the Lagrangian formulation of such equations the system 

p k = d£(x,x)/dx k 



(1.3) 

p k = -d£(x,x)/dx k 
The Hamiltonian formulation of these equations: 

x k = dH(x,p)/dp k 

(1.4) 

p k = -dH(x,p)/dx k 

The well-known Legendre transform gives a relation between the Lagrange and the Hamil- 
ton functions: 

H(x,p) = ^x k Pk ~ C(x,x) (1.5) 

k 

A discrete time analog of the Euler-Lagrange equations (|1.2|) : 

d 

(A(x, x) + A(x, x)) = (1.6) 



dx k 

We will call the Lagrangian formulation of such equations the system 

' p k = -dA(x,x)/dx k 

< 

k p k = dA(x,x)/dx k 



;i-7) 



The Hamiltonian formulation in the discrete time case is not well-defined. However, for 
the systems related to the relativistic Toda hierarchy, we will find nice approximations of 
the Hamiltonian equations of motion flO). 



2 Simplest flow of the Toda hierarchy 
and its bi— Hamiltonian structure 

The simplest flow of the Toda hierarchy (hereafter denoted by TL) is: 

(ifc = a k (b k+ i - b k ), b k = a k -a k - 1 (2.1) 

Its discretization (called hereafter dTL): 

Pk+i r , . 7 / a k a k -i\ . , 

a k = a k ——, b k = b k + h — - - — (2.2) 

Pk \Pfc Pk-\ 



where (3 k = /3k(o>, b) are defined as a unique set of functions satisfying the relations 

Pk = 1 + hbk - = l + hb k + 0(h 2 ) (2.3) 

Pfc-i 

Both the TL and the dTL are bi-Hamiltonian. The first ("linear") invariant Poisson 
bracket: 

{a k , h}i = -{a k , & fc +i}i = a k (2.4) 
The second ("quadratic") one: 



{a k ,a k+ i} 2 — —a k+ ia k , {b k ,b k+ i} 2 — —a k 
{o-k,b k } 2 = a k b k , {d k ,b k+1 } 2 = —a k b k+1 

The Hamilton functions generating the flow TL in these brackets: 



(2.5) 



Hi(a, &) = ^E^ + E a *> H ^ h ) = Y. h * ( 2 - 6 ) 

k k k 

3 Simplest flows of the relativistic Toda hierarchy 
and their bi— Hamiltonian structure 

The first flow of the relativistic Toda hierarchy (denoted hereafter RTL+): 

d k = d k (c k - Cfc_i), c k = c k (d k+1 + Cfe+i - d k - c k -i) (3.1) 
The second flow of the relativistic Toda hierarchy (denoted hereafter RTL— ): 

d\ = d k ( k 1 ] , c k = c k ( J (3-2) 

\d k d k+ i d k -\d k ) \d k d k+ i J 

An integrable discretization of the flow RTL+ (denoted hereafter dRTL+): 

t i a k+i - hd k+1 „ ofc+i + hc k+ i 

d k = d k — — , c k = c k — (3.3) 

Ofe — iia k a k + nc k 

where = a k (c, d) is defined as a unique set of functions satisfying the relations 

a k = 1 + hd k + ^zl = i + h(d k + c fc _i) + 0(h 2 ) (3.4) 

flfc-l 

An integrable discretization of the flow RTL- (denoted hereafter dRTL— ): 



where Df. = f fe(c, d) is defined as a unique set of functions satisfying the relations 

*k=j r^r. = j + 0(A) (3.6) 

An integrable difference equation (hereafter called explicit dRTL): 

J , dk + Ck ~ dk+l + Ck+l , _ x 

<4 = Gf fc _i ■ , c fc = c fc — - — ■ (3.7) 

This map is not close to the identity and therefore cannot serve as a discretization of some 
flow of the RTL hierarchy. However, it becomes related to the flow TL, if 

4 « 1 + hb k , c k « /i 2 a fc (3.8) 

(see Sect. 7 for details). 

The flows RTL± and the maps dRTL± and explicit dRTL are bi-Hamiltonian. The 
first (linear) invariant Poisson bracket: 

{ck,dk+i}i = —Ck, {ck,dk}i = Ck, {dk,dk+i}i = Ck (3.9) 

(only the non-vanishing brackets are written down), and Hamilton functions generating 
the flows ( ft.ip , (|3.2|) in this bracket: 

H{ + \c,d) = W(dk + c k ^) 2 + J2(d k + c k ^)c k (3.10) 
z k k 

H { f\c,d) = -£log(4) (3.11) 

k 

The second (quadratic) invariant Poisson bracket: 

{c/c, Cfc + i} 2 = — CfeCfe+i, {c k , d k +i}2 = —c k dk+i, {ck, d k }2 = c k dk (3-12) 
the corresponding Hamilton functions: 

H { 2 +) (c,d) = £(4 + c fc ) (3.13) 

k 

Ht\c,d) = (3-14) 

4 Systems TL and dTL: 

Parametrization of the linear bracket 

Canonical parametrization of the bracket {•, -}i: 

a k = e Xk + i-x k ^ b k = Pk (4.1) 
The resulting system is the conventional Toda lattice. 



4 



4.1 The flow TL 

Hamilton function: 



Hamiltonian equations of motion: 

Xk = Pk 



Pk = e 



Xk+l—Xk _ p Xk—Xk-l 



Newtonian equations of motion: 



x k+l~ x k _ p Xk ~ Xk-1 



x k = e 

4.2 The map dTL 

Lagrangian form of equations of motion: 

hp k = e Xk ~ Xk _ i + h 2 e x k -x k -i 

hp k = e Xk ~ Xk - 1 + h 2 e Xk + l ~ Xk 
Newtonian form of equations of motion: 



h 2 e 



Xk+l^k X k —Xk-l 

— e 



Xk X k X k X k 

e — e 

4.3 Lax representations and r— matrix structure 

The both systems above allow a Lax representation with the Lax matrix 

/ Pk + X e Xk \ 

Lk = 

\ - e ~ Xk / 
This matrix satisfies the r-matrix ansatz 

P 



{L k (X) ? Lji/M)} 



A — fj, 



>kj 



with 



/ 1 \ 

10 

10 

V 1 / 



which assures the complete integrabilty. 
The Lax representation for the flow TL: 



with 



L k — M k+ iL k — L k M k 

( -X -e x * \ 
M k = 

\ e' x k-i / 
The Lax representation for the map dTL: 



with 



L k V k — V k+1 L k 



( i-hX-h^-Xk-i - he x k \ 



V k 



V 



he-Xk-i 



5 Systems TL and dTL: 

Parametrization of the quadratic bracket 

Canonical parametrization of the bracket {•, 

a k = eZfc+i-Zfc+Pfc, bk = e Pk + e x k -x k . K 

The resulting system is the modified Toda lattice. 
5.1 The flow TL 

Hamilton function: 



H(x,p) = H 2 (a,b) =Y, ePk + ^2e x k~ x k-i 

k k 



6 



Hamiltonian equations of motion: 



Pk 



X k = 6 



= X k — g^fc %k— 1 



Newtonian equations of motion: 

x k = x k (e x k+\ — %k — e Xk — Xk-i^ 

5.2 The map dTL 

Lagrangian form of equations of motion. 

< 

he p k = (^ e x k -x k _ ^ ^ + he x k+1 -x k ^ 
Newtonian form of equations of motion: 

x k +\—x k 

_ 1 ) I i + nc 

x k x k 



(eX k -x k _ 1 ^ \l + he~ 





l + he X k~ x k-l 

5.3 Lax representations and r— matrix structure 

Lax representations for the both systems above include the Lax matrix 

/ \ e Pk _ A -i e x k 



V 



_ e -x k 



A 



satisfying the r-matrix ansatz 

{L k (X) f Lj(n)} — [r(\, /i) , L k (\) ® L k (/i) 

with 



/ 1 X 2 +fi 2 











\ 



The Lax representation for the flow TL: Q4.10|) with 

V \~ l e~ Xk - 1 



The Lax representation for the map dTL: ( 4.12 ) with 

/ I + h\~ 2 + he x k-*k-i -h\- 1 e x k\ 



V 



hX -i e -Xk-i 



} 



6 Systems TL and dTL: 

Parametrization of the mixed bracket 

For an arbitrary real e the linear combination {•, -}i+e{-, -}2 allows the following 
parametrization (a 1-parameter deformation of (|4.1| )): 



a k = e 



x k+ i-x k + ep k 



e epk _ i 
b k = h ee Xk Xk ~ 1 



6.1 The flow TL 

Hamilton function: 

H(x,p) = e-'Hzia, b) - e 1 ^Pk = £ ^ 1 \ ~ ~ + E e Xk ~ Xk ^ 

k k 6 k 

Hamiltonian equations of motion: 

' x k = (e^-l)/e 

— g-^fc+l £ k ^k—l 

Newtonian equations of motion: 

x k — (1 + ex k ) (e Xk + 1 ~ Xk — e Xk ~ Xk ~ 1 
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6.2 The map dTL 

Lagrangian form of equations of motion: 

he ep k = ^ e i^x k -x k _ 1 y j + f^j^ 1 + hee x k -x k ^ 

he^Vk = ^ (eX k -x k _ 1 ^ +h s j(i + hee x k+1 -x k ^ 

Newtonian form of equations of motion: 

e ( e x- k -x k _^ +h (l + hee~ k+l ~ Xk 



(6.5) 



x k x k 



e ^x k ±k_^ + h ( 1 + hee x k -x k _ 1 



(6.6) 



A special case e = h. If the parameter e coincides with the (small) stepsize h, then the 
previous discrete time system simplifies and serves as a discretization of the Toda lattice. 
The corresponding Lagrangian equations of motion: 



' e hp k = e x k -x k ^ + h 2 e x k -x k ^ 

< 

e hp k = e x k -x k ^ + h 2 e x k+1 -x k ^ 
Newtonian form of equations of motion: 

( -I , u2 X k+l^k\ 

-2x k +x k _{ 1 + he ~ ) 



(6.7) 



x k - 



l + h 2 e X k -X k -! 



(6.1 



6.3 Lax representations and r— matrix structure 

The both systems above have Lax representations with the Lax matrix 

/ \ e m _ A -i ee x k \ 



V 



-ee Xk 



(6.9) 



Ae 2 J 



which satisfies the following r-matrix ansatz: 

{L k (\) ® Lj(/j,)} = e\r(\,fi) , L k (\) ® L k (n) 



(6.10) 
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with the r-matrix ( |5.9| ). 

The Lax representation for the flow TL: ( f4.1(J| ) with 

/ A -2 — 1 \ 



V 



J 



The Lax representation for the map dTL: ( 4.12 ) with 



I i_^ + ^ A - 2 + /i(e-/i)e^-^-i -h\- 1 e x k\ 



V, 



e e 



V 



hx -l e -x k -l 



1 / 



(6.11) 



(6.12) 



Special case e = h. In this case also the Lax representation of the map dTL simplifies 
significantly: it reads (fO"3) with 



/ \ e hp k _ A -i he x k \ 

y -he~ Xk \h 2 j 
A" 2 -h\- l e x * 



(6.13) 



(6.14) 



7 Explicit dRTL: different parametrizations 
7.1 Hirota's discretization of the Toda lattice 

This system corresponds to the following parametrization of the (c, d) variables of the 
relativistic Toda hierarchy: 

Ck = h 2 e Xk +i~ Xk , d k = l + h Pk -h 2 e x k- x k-i (7.1) 

This results in the Poisson bracket h{-, 

The equations of motion of the explicit dRTL map in this parametrization may be put 
in the following Lagrangian form: 

hp k = e Xk ~ Xk — 1 — h 2 e Xk+1 ~ Xk + h 2 e Xk ~ Xk ~ 1 

(7.2) 

hpk = e Xk ~~ Xk - 1 
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It follows a nice " Hamiltonian" form of equations of motion: 



e Xk Xk — 1 = hp k 

p k —p k = h^k+l—Xk _ fo e Xk— Xk-l 



(7.3) 



The Newtonian form of equations of motion: 



%k %k %k %k 



h 2 e 



x k+l~ x k x k ~ x k—\ 



(7.4) 



The Lax representation: (|4.12 ) with 

/ p k + A e Xk \ 

Lk = 

\ -(1 + hp k )e~ Xk -h J 



I 1 - hX - h 2 e Xk ~ Xk - 1 -he Xk 



V, 



V 



h e - Xk -l 



(7.5) 



(7.6) 



1 



Note that the matrix V k depends only on the x/s, and not on their discrete time updates 
Xj or on the momenta pj. This will be the common feature of all the results in this Section. 



The Lax matrix ( |7.5| ) is a one-parameter deformation of the standard Toda Lax matrix 
Q4.7D , but still satisfies the r-matrix ansatz (|4.8|) 

7.2 Standard discretization of the Toda lattice 

This system corresponds to the following parametrization of the variables (c, d) of the 
relativistic Toda hierarchy: 

Ck = h 2 e x k+ i-x k + hp kj d k = e h P k (7.7) 

This results in the Poisson bracket h{-, -}2- The equations of motion of the explicit dRTL 
map may be presented in this parametrization in the following Lagrange form: 



JiPk 



,hpk 



(l + h 2 e X k -X k ^ 



■y^k *^k 



1 + h 2 e Xk+1 ~ Xk 



(7.1 
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This may be presented in a nice " Hamiltonian" form: 

x k - x k = hp k 



(l + h 



The Newtonian form of equations of motion: 

x k -2x k + x k _ (l + h 2 e x k+i~ x k 



(l + h 2 e Xk ~ Xk - 1 

The Lax representation: (|4.12 ) with 

( \ e hp k _ A -i he x k \ 

L k = 



( 



V k 



_ he -x k +h Pk o ) 
A- 2 -h\- l e Xk 



(7.9) 



(7.10) 



(7.11) 



(7.12) 



V hX^e-^-i 1 

The Lax matrix ( |7.11| ) is also a 1 -parameter deformation of the standard Toda Lax matrix 
( |4.7| ), but satisfies a different r-matrix ansatz: 

{L k (X) f Lj(fA)} = h[r(X,n) , L k {\) ® L k (ji)]6 kj (7.13) 

with the r-matrix 



7.3 Explicit discretization of the modified Toda lattice 

This system corresponds to the following parametrization of the variables (c, d) of the 
relativistic Toda hierarchy: 

Ck = h 2 e Xk + 1 ~ Xk+ P k , 4 = 1 + heP k + he Xk ~ x ^ (7.14) 

This results in the Poisson bracket {•, — {•, The equations of motion of the explicit 
dRTL map may be presented in this parametrization in the following Lagrange form: 



heP k 



. (l + he Xk ~ Xk -^) 
) (l + he Xk +i- Xk ) 



(7.15) 



^ e x k -x k _ A 
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This may be presented in a nice " Hamiltonian" form: 

f e x k -x k _ 1 = he p k 



(l + ht 



e Pk~Pk 

(l + he Xk ~ Xk ~ 1 
The Newtonian form of equations of motion: 



1 + he Xk Xk ~ 1 



The Lax representation: (|4.12 ) with 

/ \ e Pk _ a- 1 

Lk = 

V 



V, 



l + heP k )e~ Xk A 
/ 1 + ii\- 2 + he Xk ~ Xk -i //A ; r 

V 



/iA _1 e -a; *-i 1 
The Lax matrix (738 ) is a 1-parameter deformation of the Lax matrix ( |5.7| ), but 



the same r-matrix ansatz (5.8 



8 Systems RTL± and dRTLi: 

Parametrization of the quadratic bracket 

Canonical parametrization of the bracket {•, 



di 



oPk 



Ck 



g 2 e X k +l—Xk+Pk 



The resulting system is the relativistic Toda lattice. 
8.1 The flow RTL+ 

Hamilton function: 

H{x,p) = H ( 2 + \c,d) = Y / eP k (l + g 2 e Xk +i- Xk 
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Hamiltonian equations of motion: 

' ± k = eP k (l + g 2 e x k+i- x k^j 

< 

^ p k = g2 e X k+ l-X k +Pk - g 2 e Xk-X k -l+Pk-l 

Lagrangian equations of motion: 

e Pk 



x k 



(l + g 2 e X k+1 -X k ^ 



g2^X k +\ X k g2^X k X k _\ 

Pk = X k — ^ ^ x - X k -1 



(l + g 2 e X k+1 -X k ^ (\ + g 2 e X k -X k _i^ 

Newtonian equations of motion: 

g2{^Xk+\ X k g2^X k X k —\ 

X k X k J r \X k — - — X k X k —\ 



(l +g 2 e X k+1 -X k ^ 



2 p X k -X k -i 



g*e 



) 



8.2 The map dRTL+ 

"Hamiltonian" form of equations of motion: 

e x k — x k fogPk^i _|_ g2 e x k+ \— x k ^j 



1 + hg 2 e Xk + 1 ~ Xk 



+Pk\ 



e Pk - Pk 

fl + hg^k-Xk-l+Pk-l 
Lagrangian form of equations of motion: 

f e Xk ~X k _A fl + g 2 e X k 

h(l'>- 



I _|_ g2 e x k— Xk-1 j ( [ . gl ( •''/■• • I 

g-^fc Xk ~y 



heP k = 



I _|_ g2 e Xk + l~Xk 
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Newtonian form of equations of motion: 



Xk—Xk _ I j ( 1 + g 2 e Xk + 1 ~ Xk ) ( 1 + g 2 e Xk ~ Xk ~ 1 



Xk Xk ^ 

e ~ — 1 



„ Xk+l—Xk\ I o Xk—Xk-l 

1 + c/ 2 e~ 1 + g 2 e 



8.3 The flow RTL— 

Hamilton function: 

H(x,p) = H { 2 ~\c,d) = J2e~P k (l + g 2 e Xk ~ Xk -^ 

k 

Hamiltonian equations of motion: 

x k = -e~P k (l+g 



2gXfc Xk— 1 



p k — g 2 e Xk+l ~ Xk~ Pk+1 _ g 2 £ X k — Xk-1 — Pk 



Lagrangian equations of motion: 



oPk 



(}+g 



I^Xk Xk—1 



Pk — —Xk+l 



X k 



g 2 ff^k-\-l Xk 



+ X k 



g 2 (^Xk Xk—\ 



(l + g 2 e Xk+1 ~ Xk ^j (l + g 2 e Xk ~ Xk ~ 1 
Newtonian equations of motion - the same ( ]8.5p as for the RTL+ (!). 

8.4 The map dRTL— 

"Hamiltonian" form of equations of motion: 

oXk -x k _i = -he~P k (l + g 2 e Xk ~ x k-i 



oPk - Pk 



I — hg 2 e Xk—Xk-l—Pk 
I — f l g2 e Xk+\—Xk—pk+l 
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Lagrangian form of equations of motion: 

h ( 1+g 2 e X k -X k - 1 

ppk = -- 



oPk 



£>X k X k \ 

h(l + g 2 e Xk + 1 ~ Xk ) ( 1 + g 2 e Xk ~ Xk - 1 



(8.13) 



Newtonian form of equations of motion - ( |S.8|) , the same as for the dRTL+ (!). 

8.5 Lax representations and r— matrix structure 

Lax representations for the four systems of this Section is given in terms one and the same 
Lax matrix: 



\ e Pk _ A -i e x k 



-q 2 e - x k+Pk o 



(8.14) 



It satisfies the r-matrix ansatz ( |5.8j ) with the r-matrix ( p.9| 
Lax representation for the flow RTL+: (^4.10|) with 



/ ^-2 _|_ q "2 e Xk-Xk-l+Pk-l -\-i e Xk \ 



M, 



Vk 



\ 



V \-i g 2 e -x k -i+Pk-i o J 
Lax representation for the map dRTL+: ( [4.12 ) with 

j i + h \-2 + hg 2 e Xk-x k -i+Pk-i -h\- l e Xk \ 

hX~ l g 2 e- Xk -i+Pk-i i 

Lax representation for the flow RTL— : (|4.10 ) with 

/ o -\ e x k -Pk 

M k = 

V \g 2 e- x *-i X 2 + g 2 e Xk ~ Xk ^ ~P k 
Lax representation for the map dRTL— : 

W k+1 L k = L k W k 



.15) 



.16) 



(8.17) 



Hi 
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with the same Lax matrix ( |3.14|) and 

/ 1 



hXe X k~Pk 



i -hXg 2 e Xk ~ 1 1 - hX 2 - hg 2 e 



9 Systems RTL± and dRTLi: 

Parametrization of the linear bracket 

Canonical parametrization of the bracket {•, 

d k = p k — e Xk ~ Ck = e Xk + 1 ~ Xk 

9.1 The flow RTL+ 

Hamilton function: 



H(x, p) = H[ +) (c,d) = 5Erf + E P keXk+1 

1 k k 



Hamiltonian equations of motion: 

x k = p k + e Xk+1 ~ Xk 

p k = p k e Xk +i~ Xk — p k _ 1 e Xk ~ Xk ~ 1 
Lagrangian equations of motion: 
p k = x k — e Xk + 1 ~ Xk 

— (^x k Q X k+l x k^j p^k+1 X k {^X k 1 P Xk X k—l^j Q X k x k 

Newtonian equations of motion: 



Xk = [Xk+i - e 



_ p x k+l ~ x k 1 P x k+l~ x k 



Xk-i - e 



p X k x k—l ) p X k ' 



17 



9.2 The map dRTL+ 

"Hamiltonian" form of equations of motion: 

e x k -x k _i = j l p k + 



)i e x k+l ~ x k 

1 - he x k+i~^k 



(9.6) 



Pk - Pk = hp k e Xk +^ Xk - hp k _ 1 e Xk Xk ~ 1 
Lagrangian form of equations of motion: 

1 



hpk = e Xk Xk 



hp k = e Xk Xk — 



1 - he Xk ~ Xk -± 
1 



1 - he Xk +i~ Xk j 



+ he Xk Xk ~ l — he Xk+1 Xk 



(9.7) 



Newtonian form of equations of motion: 



x k x k x k x k ^Cfc+i x k he 
e — e ~ = he 



X k+ i^C k 



I -he 



Xk+l^Ck 



- he Xk Xk ~ l + 



he x k -x k ^ 

1 - he Xk ~ Xk ~i 
(9.8) 



9.3 The flow RTL— 

Hamilton function: 



H(x,p) = Ht\c,d) = -E lo g (Pk - e^-^-i) 

k 

Hamiltonian equations of motion: 



(9.9) 



Pk 



{j>k — e Xk Xk -^ 
e x k +i — x k 



D X k X k —\ 



(9.10) 



(pk+i — e Xk + l Xk ^j (p k — e Xk Xk ~^j 
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Lagrangian equations of motion: 



p k = -— + e X k~ x k-l 
x k 

k p k = -x k+1 e Xk + 1 ~ Xk +x k e Xk ~ Xk - 1 
Newtonian form of equations of motion: 

x k = -x 2 k {x k+1 e Xk + 1 ~ Xk - Xk-ie Xk ~ Xk - 1 ) 

9.4 The map dRTL— 

"Hamiltonian" form of equations of motion: 

^>X k Xk ~y — _ 



h 



Pk ~Pk 



Pk — e Xk Xk ~ l 

flQ-Ek+l Xk h,C Xk "^k— 1 



Pk+i - e Xk+1 Xk J [ Pi ' 



Xk Xk—1 



Lagrangian form of equations of motion: 

h 



Pk 



Pk 



gXfc Xk J 

h 

?x k x k 



_|_ gXk Xk—1 



_|_ e x k+l x k _ e X k+ i Xk _|_ e Xk X k -1 



Newtonian form of equations of motion: 



h 



h 



(9.11) 



(9.12) 



(9.13) 



(9.14) 



,X k Xk ~y 



Xk Xk 
e ~ — 1 



x k+l~Xk Xk+l—Xk Xk~Xk-l Xk 

e — e~ — e + e 



— Xk-l 



(9.15) 

9.5 Lax representations and r— matrix structure 

All four systems considered in this Section allow Lax representations with the Lax matrix 

f p k + \ e Xk \ 
L k = (9.16) 
V ~p k e~ Xk -1 / 
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which satisfies the r-matrix ansatz (4.8). 

Lax representation for the flow RTL+: (4.1U) with 



-A 



pXk 



Pk-\e 



-Xk-l 



(9.17) 



Lax representation for the map dRTL+: fl4.12| ) with 

( \-h\-h 2 Vk-\e Xk ~* k - x -he Xk \ 



hp k _ 1 e~ Xk ~ 1 

Lax representation for the flow RTL— : ( |4.10|) with 

Pk 



(9.is; 



Mk 



A 



-i 



oXk 



-Pke 



(pk — e Xk Xfc_1 
Lax representation the map dRTL— : ( 8.18|) with 

h 1 



Xk— 1 pX k X k — 1 



(9.19) 




This last Lax representation may be also presented in the form ( |4.12|) with 
being more nice than W k (which is non-tipical for the map dRTL—): 



V k = I + 



hX 



Pk e Xk 
p k - e Xk-x k -i) { e -x k . x _ e x k -x k _ 1 



(9.20) 
matrix V k 

(9.21) 



10 Systems RTL± and dRTLi: 

The first mixed parametrization 

We now consider the following canonical parametrization of the bracket {-, — 5{-, •}].: 



di 



aPk 



i p x k -x k -i 



Ck = g 2 e Xh+l-X k +p k 



10.1) 



This is, obviously, a 1-parameter deformation of (jTT). The results of this section in the 
limit 5-^0 imply the results of the Section 8. 
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10.1 The flow RTL+ 

Hamilton function: 

H(x,p) =H { 2 + \c,d) = £eP*(l + 2 e a *+i- a; *) + 5g 2 ^ Xk ~ Xk - 1 

k k 

Hamiltonian equations of motion: 

± k = eP k (l + g 2 e X k+l~ x k\ 

Lagrangian equations of motion: 



(10.2) 



(10.3) 



X k 



Pk = x k 



(l + g 2 e Xk+1 ~ Xk 

g 2 k+1 Xk 



(1+ 



2 p X k+ i — X k 



g*e 



- X k -l 



g 2 gX k Xk—1 

(l + g 2 e Xk ~ Xk - 1 



§g 2 (g^k-{-l X k ^X k X k —i 

(10.4) 



Newtonian equations of motion: 



x k 



Xk+lXk 



g2^X k -\-\ X k 



(l _|_ g2 e X k + l—X k 

+5g 2 x k (e Xk + 1 ~ Xk -e x k~x k -\ 



XkXk—l 



g 2 (^Xk X k —\ 
(l + g 2 e Xk ~ Xk ~ 1 



10.5) 



It is interesting to remark that the right-hand side of this equation coincides with the 
right-hand side of flS.5| ) plus an additive perturbation which is exactly the right-hand side 
of the modified Toda lattice 05.41). 



10.2 The map dRTL+ 

"Hamiltonian" form of equations of motion: 

heP k (l + g 2 e Xk + 1 ~ Xk 



(>Xk X k 



e Pk~Pk 



l + h5g 2 e Xk + 1 ~ Xk 
l + hg 2 (eP k +§y x k+i-Xk 



(10.6) 



(l + hg 2 (ePk-i + 8)< 



•Xk X k —\ 
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Lagrangian form of equations of motion: 

f e X k -X k _l\ (l + g 2 e x k-X k -l 

hePk = 



I _|_ g 2 e Xk ~ X k-l ) ( 1 4- g1 e Xk+l - ' 
gXk X k ^ 



1 + h5g 2 e x k~ x k-i 



(10.7) 



l + h5g 2 e x ^- x k 

I _|_ g2 e Xk + l—Xk \ V / 



Newtonian form of equations of motion: 

Xk X k 



1 



1 + g 2 e Xk+l-X k g 2 e X k -X k „ 1 ) ( 1 + ^2 e ^+l^A 



e ~ — 1 



(10.8) 



Interesting enough, the right-hand side of the last equation is the product of the right- 
hand sides of (§3) and 



10.3 The flow RTL— 

Hamilton function: 

H(x,p) = -5- l H[-\ Cl d)-5- l Y,Pk = r x ^log(4) -S~ iy £Pk 

k k k 

where, recall, 

ePk + 0(5) 



4 = e Pk + 5(1 + ^2 e x fc -x fc _i 
Hamiltonian equations of motion: 

Xk = -d^il + g 2 e X k- x k-l 



(10.9) 
(10.10) 

(10.11) 



Pfc = d k l 1 g 2 e Xk + 1 Xk - d k 1 g 2 e Xk Xk ~ l 
Lagrangian equations of motion: 



oPk 



(l +g 2 e X k -X k ^ ^ 



Vk — ~ X k +1 



x k 

g 2 qX k+1 X k 
2 p X k+ i-X k 



5x k ) 



+ x k 



g 2 (f^k X k —\ 



(10.12) 



g*e 



g2^X k X k —\ 
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Newtonian equations of motion: 



x k = (1 + 8x k ) x k+ iik 



g 2 e %k + l ~X k 



2 p X k+ i-X k 



9 e 



- X k X k -i 



g 2 p^k X k — .1 



2 p X k -Xk-l 



9 e 



This is a multiplicative perturbation of 



(10.13) 



10.4 The map dRTL— 

" Hamiltonian" equations of motion: 

' e Xk~Xk _ 1 = -hD^(l + g 2 e X k-^k-l 



oPk ~Pk 



l-hD^g^k-Xk-i 



1 - hD^ ig 2 e x k+i 



-x k 



where 

D k = ePk + 5(l + g 2 e x k~ x k-i^ = e Pk + 0(5) 
Lagrangian form of equations of motion: 



hePk 



5^k-x k _ +h 



-,X k X k j 



X ~\- n 2 e x k~x k _i 



heP k 



-,X k X k j 



+ g 2 e Xk+ ^ Xk 



\ _|_ g 2 Q X k X k —\ 



Newtonian form of equations of motion: 



l\U(e Xk ~ k -l) + h} (l+g 2 e x *+i- x ^(l + g 2 e 



;io.i4) 



;io.i5) 



(10.16) 



Xk X k —\ 



x k x k 



-l)(6(e x k- x k-l)+h\ (l+g 



i2f 3+i + /y2r r,-.r,_j 



;io.i7) 
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Special case 5 = h. If the parameter 6 is equal to the (small) stepsize h 
dRTL— simplifies, providing discretization of the relativistic Toda lattice 
from ( |8.8|) . In this case the " Hamiltonian" equations of motion read: 



1 



1 + he~P k [l + g 2 e Xk ~ 

( 1 + hg 2 (ePk+i + h) ~ -x k \ 
e p k -Pk = V ^ '_ / 

(l+hg 2 (fk+hy 1 e c >>- x >>-^ 
The Lagrangian equations of motion read: 

h(\ ■ ,/V/."-' 

oVk = _ . 



1-e" 



oPk 



h(l + g 2 e Xk ~ Xk - 1 ) ( 1 + g 2 e Xk + 1 ~ Xk 



■x k 



The Newtonian equations of motion: 
1 — e 



1 + g 2 e Xk+i-x k j ( l + ^e 



Xk—1 



1-e 



-%k~\~ Xk 



J (l + </V 



the map 
different 



'10.18) 



;i0.19) 



(10.20) 



The system resembles very much the previous discretization of the relativistic Toda lattice 
~M, however the relation between the two is far from trivial. 



10.5 Lax representations and r— matrix structure 

The Lax representations for the four systems considered in this Section may be given in 
terms of the following Lax matrix: 



( XePk-X- 1 e Xk \ 
_ g 2f e p k + S ) e -Xk \5g 2 



'10.21) 
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which satisfies the r-matriz ansatz (|5.8| ) with the r-matrix ( |5.9|) . 
The Lax representation for the flow RTL+: ( 4.10|) with 



M k = 







Lax representation for the map dRTL+: ( |4.12 ) with 

/ 1 + hX~ 2 + hg 2 (fk-\ + §\ e Xk-Xk-i -hX-^k \ 

V k = 

v h\- l g 2 (ePk-i + $y-Xk-i 
Lax representaion for the flow RTL— : (|4.10| ) with 

^ Xg 2 (eP k + 5)e~ x >*- 1 g 2 e x k-x k -i 
Lax representation for the map dRTL— : (|8.18 ) with 

/ \ 2 (ePk + 5) Xe Xk \ 



(10.22; 



(10.23) 



d' 1 
k 1 + 5X 2 



-Xe Xk \ 
) 



(10.24) 



W k = I + 



l + (5-h)X 2 



[ -Xg 2 (ePk + 5)e- Xk -i -g^^k-x^ 



(10.25) 



Special case 5 = h. In this case, for the map dRTL—, the dependence of the matrix 
W k on A simplifies, because the denominator 1 + (5 — h)X 2 becomes equal to 1. 

11 Systems RTL± and dRTLi: 

The second mixed parametrization 



We consider the following canonical parametrization of the variables (c,d): 

e ep k _ 1 



di 



^X k X k —\ 



,x k+ i-x k + ep k 



11.1) 



By small e it serves as a deformation of (|9.1| ). The corresponding Poisson bracket is 
{•, -}i + e{-, -}2. the results of this Section allow to recover that of the section 9 in the 
limit e — > (in order to reproduce the Lax representations of the Sect. 9, the spectral 
parameter A of this Section has to be replaced by 1 + eA/2 + 0(e 2 ) before perforimg the 
limit e — > 0). 
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11.1 The flow RTL+ 

Hamilton function: 

H(x, P ) = e' l H^ (c, d) -e- 1 J2Pk = J2 ~ + E " ~ 

k k e k e 

Hamiltonian equations of motion: 

(e € P* - l) 



x k 



Pk 



+ e x k+ i-x k + ep k 



UPk-i) ( e ePk-i-i) 
A LpXk+i-Xk _ V )_ 



Q X k %k—l 



Lagrangian equations of motion: 
(1 + ex k ) 



^Pk 



Pk 



+ l X k 



x k — e 



g^fe+l x k 



( 



x k -i — e Xk Xk ~ 1 



x k x k 



(l + ee x k+ i-x k ^ (i + ee x k -x k -^ 

The corresponding Newtonian equations of motion read: 



-0 



'x k = (1 + ear*.) 



(l + ee Xfc + 1_Xfc ) 



(l + €e x k -x k _^ 



?x k 



11.2 The map dRTL+ 

"Hamiltonian" equations of motion: 



( e €Pk-i) h e x k +i-x k + ep k 

e x k -x k _ x = h \ ; + 



1 - ^fc+i-aJfcj 



e ep k -ep k = 



(l + h(e e Pk -ly^k+i-x^ 
(l + h(e € Pk-i -l) e ^-^-i) 



26 



Lagrangian equations of motion: 



I _ he Xk ~ Xk ~ 1 ) (i 4- te Xk ~ x k-i 



he e P k = e e Xk ~ Xk -l)+h 



1 + (e-h)e Xk ~ x k-i) [i + ee Xk+i-x k 

I — f le X k +l~X k 



11.7) 



1 + (e- h)e Xk + 1 ~ Xk 



Newtonian equations of motion: 

X k X k 



e e 



l) +h (i + ee x k+ i-x k ) {i-he~ k+1 ^ k ) {i + ( e - h )eCk-Xk-i 



c\e Xh ~ k -l)+h (l + ee Xk - Xk -A(l-he Xk - Xk -A(l + (e-h)e~ k+1 %k 



11.8 



Special case e = h. If the parameter e is equal to the (small) stepsize h, then the 



equations above greatly simplify, delivering another discretization of the system (9.5). In 
this case we obtain the following " Hamiltonian" equations of motion: 



?x k x k 



1 - he Xk +i~ Xk ) 



MPk-hp k 



Lagrangian equations of motion: 



1 + h (ehp k _ 1 yx k+1 -x k \ 
fl + h(e h Pk~i - l)e Xk ~ Xk ~^J 



11.9) 



Jipk = e Xk~X k {]__ foXk-Xk-AS^ 

^ Ml 



1 + he Xk Xk ~ l 



;n.io) 



,hp k _ e X k — X k M _ fo e Xk+l — x k 
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Newtonian equations of motion: 

^_ 2li+Jl _(i + ^«-^)(i-/ K ?-^) 



1 + he x k~ x k-i 1 - he Xk ~ Xk ~ l 



11.3 The flow RTL— 

Hamilton function: 

H(x, p) = H[~ ] ( c ,d) + eJ2Pk = l °g( d k) + e ^p k 

k k k 

where, recall, 

e ePk _ i 

4 = e Xk Xk ~ l = p k - e Xk Xk ~ 1 + 0(e) 

Hamiltonian equations of motion: 

' x k = -d^ (l + ee Xk ~ Xk -^ 

dzL e Xk + l ~ Xk - dZ 1 e Xk ~ Xk ~ l 



{ Pk = %+i 



Lagrangian equations of motion: 



£Pk = (i + ee x k 



-x k -,\ ( Xk ~ £ ) 

Xk 



(yX k -\-\ X k ^X k X k _i 

Pk = —x k+ \ — — + x k 



(l + ee Xk + 1 ~ Xk ) (l + ee 2 * -3 *- 1 ) 



Newtonian equations of motion: 
x k = —x k (x k — e) x k+ i 



X k — i 



(l + ee ^+i-^) 



(l + ee a*-a*-i) 



11.4 The map dRTL— 

"Hamiltonian" equations of motion: 

' e x k -x k _ 1 ^_ hD -ifi + ee x k -x k - 1 ^ 



e ep k -ep k = 



1 - heDZ x e Xk ~ Xk ~^ 
1 - heD k l ie Xk +±- Xk 
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where 



e epk _ i 

D fe = e x k -x k ^ =pk _ e Xk~Xk-i + o( e ) 



Lagrangian form of equations of motion: 
e Xk ~ Xk — I — he 



x k 



I _|_ gg^fc - Xk-l 



e Xk ~ Xk -l-he) , ^ x [\ + ee Xk-Xk-i 

1 _)_ gg^fc+1 ~~ ^fc 



X -)- gg^fe+1 ^Tft 



Newtonian form of equations of motion: 



,x k -x k _ 1 \ t Xk 5 fc _i_/, e 



\ _|_ gg^fc+l X k \ []__)_ gg-^A; 



l + ge^+^M [i + ee x k -x k ^ 



:n.i8) 



;n.i9) 



(11.20) 



11.5 Lax representations and r— matrix structure 

The Lax representations for the four systems considered in this Section may be given in 
terms of the following Lax matrix: 

/ \ e ePk _ a- 1 ee Xk \ 

U = (11.21) 
Lep k -i) e -Xk _gA 



V 



/ 



which satisfies the r-matriz ansatz ( 6.10|) with the r-matrix (5.9). 
The Lax representation for the flow RTL+: ( |4.10j ) with 



/ X - 1 



_l_ f e ep k -! _ i^ e x k -x k ^i _\-i e x k 



V 



A 



i ((fiPk-i _ i 



-Xk-l 



Lax representation for the map dRTL+: (|4.12 ) with 

e e 



+ -A- 2 + ^ — 'll( e m-i - i) e Xk-x k -! _ h \-i e x k \ 



V 



'11.22) 



;il.23) 
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Lax representaion for the flow RTL— : ( 4.10Q with 

/ A 2 (e e ^-1 



M k 



A 2 - 1 



\ee Xk \ 



-ee 



;n.24) 



Lax representation for the map dRTL— : (|8.18 ) with 

/ A 2 (e e ^-l) Aee 2 * \ 

v -\UPk - l) e -3fc-i - e e x k-Xk-i J 



W k = I + 



1 - (1 + he)X 2 



;il.25) 



Special case e = h. In this case the Lax representation for the map dRTL+ simplifies 
significantly: it reads ( JQ2Q with 



Vk 



( Xe h Pk - A" 1 he x k \ 

e hp k _ i) e -x k _h\ j 

A- 2 -hX- 1 ^ \ 

K X- l (e h Pk-i-l) e - x k-i i 



;il.26) 



;n.27) 
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